CUBIC THUE EQUATIONS 



SHABNAM AKHTARI 



Abstract. We revisit a work by R. Okazaki and prove that for every cu- 
bic binary form F(x,y) with large enough discriminant, the Thue equation 
\F(x,y)\ = 1 has at most 7 solutions in integers x and y. 



1. Introduction 

Let F(x, y) be an irreducible binary cubic form with integral coefficients and 
negative discriminant. More than 80 years ago, Delone and Nagell established 
independently that the equation 

(1) \F(x,y)\ = l 

has at most five solutions in integers x, y. This result is proved by considering units 
in the algebraic number field Q(p), where p is the real root of F(x, 1) = 0. In their 
proofs the fact that the group of units in the ring of integers of Q(p) is generated 
by one fundamental unit is essential. 

The situation where the discriminant of F(x, y) is positive is complicated by 
the fact that the number field Q(p) (where p is any real root of F(x, 1) = 0) 
has a ring of integers generated by a pair of fundamental units. However, it is 
possible to reduce |T]) to a set of exponential equations to which a local method of 
Skolem can be applied. In this way, Ljunggren [14] and Baulin [4], solved |T]) for 
F(x,y) = x 3 — ixy 2 + y 3 of discriminant 81 and F(x,y) — x 3 + x 2 y — 2xy 2 — y 3 
of discriminant 49, respectively. In the first case there are 6 solutions and in the 
second case there are 9 solutions to (P). 

In 1929, Siegel [TB] used the theory of Pade approximation to binomial functions 
(via the hypergeometric functions), to show for F cubic of positive discriminant, 
that equation |lj has at most 18 solutions in integers x and y. Refining these 
techniques, Evertse [5] reduced this upper bound to 12. Later, Bennett [2] showed 
that if F(x, 1) has at least two distinct complex roots, then the equation F(x, y) = 1 
possesses at most 10 solutions in integers x and y. In 2003, by studying the geometry 
of numbers in the "logarithmic space" , Okazaki [15j proved that if discriminant of 
F is greater than 5.65 x 10 65 then equation has at most 7 solutions. Okazaki's 
method is essentially different from Evertse's. In this paper, we will relate some 
geometric ideas of Okazaki [15] to the method of Thue-Siegel as refined by Evertse 
[5], in conjunction with lower bounds for linear forms in logarithms of algebraic 
numbers. The following are the main results of this paper: 
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Theorem 1.1. If F(x, y) is a binary cubic form with discriminant D > 1.4 x 10 57 , 
then the equation 

\F{x,y)\ = l 

possesses at most 7 solutions in integers x and y. 

Theorem 1.2. Let F(x,y) be a cubic form with discriminant D > 9 x 10 58 . 
// F(x, y) is equivalent to a reduced form which is not monic, then the equation 
F(x, y) = 1 posses at most 6 solutions in integers x and y. 

Despite the numerical improvement, the bounds remain out of reach of comput- 
ers. The main purpose of this paper, is to look at some beautiful geometric ideas 
of Okazaki [15j from Classical Analysis point of view. 

In 1990, using the fact that the underlying number fields are the so-called "sim- 
plest cubics" , Thomas [T7| showed that the equations 

Gi, n (x, y) = x 3 + nx 2 y - (n + 'i)xy 2 + y 3 = 1 

have only the solutions (1, 0) , (0, 1) and (—1, —1) in integers, provided n > 1.365 x 
10 7 . This restriction was later removed by Mignotte [11] except for the equation 
with 

n e {-1,0,2}. 

It is known that Gi >n (x,y) = 1 has 9 solutions for n = — 1 (_!]), 6 solutions for 
n = ([14]) and 6 solutions for n = 1 ([6]) . 
Define F m {x,y) by 

F m {x, y) = x 3 - (to + l)x 2 y + mxy 2 + y 3 

for m G Z. Provided m ^ —2 , —1 or 1 , the equation F m (x,y) = 1 has the five 
distinct integral solutions (x,y) = (1,0) , (1,1) , (1,— m— 1) , (0,1) and (m, 1). 
That this list is complete was proven, independently, by Lee [7] and Mignotte and 
Tzanakis [13], for m suitably large and later, by Mignotte [12] , for to > 2. The 
cases to = and to = 1 correspond to discriminant —23 and —31, respectively. 

All known irreducible cubic forms F(x,y), for which the equation (fT]) has more 
than 5 solutions, have discriminant less than 362. 

The following conjecture is essentially due to Nagell and refined by Petho and 
Lippok. 

Conjecture If F is a binary cubic form with positive discriminant Dp, then the 
number of solutions of equation |T]) is less than 6, if Dp > 361. 

2. The Covariants of Binary Cubic Forms 

Suppose 

F = ax 3 + bx 2 y + cxy 2 + dy 3 

with discriminant 

Dp = 18abcd + b 2 c 2 - 21a 2 d 2 - Aac 3 - Ab 3 d = a 4 ]J(a, ; - a 3 ) 2 

where ai, ai and «3 are the roots of polynomial F(x, 1). Let us define, for the form 
F, an associated quadratic form, the Hessian H = Hp, and a cubic form G = Gf, 



CUBIC THUE EQUATIONS 



3 



by 

m \ = 1 ( s2fs2f f s2F Y\ 

(XlV> 4 \Sx 2 Sy 2 \SxSy) J 

= Ax 2 + Bxy + Cy 2 

and 

G(x y) - SF 6H SF SH 
' Sx Sy Sy Sx 

These forms satisfy a covariance property; i.e. 

Hfoj = Hp o 7 and Gf 01 = Gf ° 7 

for all 7 G GL 2 (Z). 

We call forms F\ and F 2 equivalent if they are equivalent under GL2(Z)-action; 
i.e. if there exist integers a± , a 2 , a 3 and 0,4 such that 

Fi(aix + a 2 y, a 3 x + a 4 y) = F 2 (x,y) 

for all x, y , where 0104 — 0203 = ±1. 

We denote by Np the number of solutions in integers x and y of the Diophantine 
equation |T]). li Fi and F 2 are equivalent, then Np 1 = Np 2 and Dp 1 = Dp 2 . 
Therefore, we can assume that F is monic (the coefficient of x 3 in F(x,y) is 1). 

For F{x, y) — ax 3 + bx 2 y + cxy 2 + dy 3 with discriminant D, it follows by routine 
calculation that 

A = b 2 - 3ac, B = bc- 9arf, C = c 2 - 3bd 

and 

B 2 - AAC = -3D. 
Further, these forms are related to F(x, y) via the identity 
(2) 4H(x,y) 3 = G(x,y) 2 + 27DF(x,y) 2 

Binary cubic form F is called reduced if the Hessian 
H(x, y) = Ax 2 + Bxy + Cy 2 

of F, satisfies 

C > A > \B\ 

It is a basic fact (see [I]) that every cubic form of positive discriminant is equivalent 
to a reduced form F(x,y) . The reader is directed to [1] (chapter III and supplement 
I) for more details on reduction. We will later use the following lemma to bound 
the discriminant D from above. 

Lemma 2.1. Let F be an irreducible binary cubic form with positive discriminant 
D and Hessian H. For all integer solutions (xi,yi) to equation F(x,y) — 1, except 
possibly one solution, we have H(xi,yi) > i %/3-D. 

Proof. If Fi is an equivalent reduced form to F and 

Fi(oia; + a 2 y, a 3 x + a 4 y) = F(x, y), 

then 

Hi{a\x + a 2 y, a 3 x + a±y) = H(x,y), 
where H and H\ are the Hessians of F and Ft respectively. This means the set of 
values of the Hessian at solutions is fixed under GL2(Z)-action. So we may assume 
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that F is reduced. Now following the proof of lemma 5.1. of [2], we suppose (x,y) 
is a solution to F(x,y) = 1 with y ^ 0. If \y\ < \x\, then, since A > \B\ and 
B 2 - A AC = -3D, we have that 

(x, y) = Ax 2 + Bxy + Cy 2 > Cy 2 > C > \V3D. 

If, on the other hand, \y\ > \x\ + 1, then 

H(x,y)> (C-\B\)y 2 + \B\\y\ + Ax 2 . 
Since this is an increasing function of \y\ and y ^ 0, we have 

H(x, y)>C + Ax 2 >C> ^V3D. 

Therefore, if H(x, y) < ^ V3D, then y ~ and so x = ±1 accordingly . □ 

Remark. The above proof shows the only possibility for the Hessian H (x, y) to 
assume a value less than ^^/3D, at a pair of solutions (x, y), is when the equivalent 
reduced form is monic. This is because (1, 0) is a solution to ([T|) if and only if F is 
monic. 

3. Some Functions In The Number Field Q(V-3D) 



Let \J — 3D be a fixed choice of the square-root of —31?. we will work in the 



number field M = Q(\/—3D). It is well-known that if F has positive discriminant 
then H is positive definite. By (|2|), we may write 

H(x,y) 3 = U(x,y)V(x,y) 

where 

_ G(x, y) + 3V^3DF(x, y) 

U \ x i V) — ,-, i 



G{x,y)-3^r3DF( Xl y) 
y{x,y) = . 

Then U and V are cubic forms with coefficients belonging to M such that corre- 
sponding coefficients of U and V are complex conjugates. Since F must be also 
irreducible over M, U and V do not have factors in common. It follows that U(x, y) 
and V(x, y) are cubes of linear forms over M, say £(x, y) and T](x, y). 

Note that £(x, y)ij(x, y) must be a quadratic form which is cube root of H(x, y) 3 
and for which the coefficient of x 3 is a positive real number. Hence we have 

£{x, y) 3 - v(x, y) 3 = 3^3DF{x, y), 

(3) C(x,y) 3 +v(x,y) 3 = G(x,y), 

£{x,y)ri(x,y) = H(x,y). 



and 



£(^1 v) , v(x, y) 

and ') y ' e M. 



£(1,0) ' »/(l,0) 

The reason for the last identity is that for any pair of rational integers xq , yo, 

£{x ,yo) and ri(x ,y ) 

are complex conjugates and the discriminant of H is —3D. 
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We call a pair of forms £ and rj satisfying the above properties a pair of resolvent 
forms. Note that there are exactly three pairs of resolvent forms, given by 

where w is a primitive cube root of unity. 

We say that a pair of rational integers (x, y) is related to a pair of resolvent forms 

if 



(4) 



1 



Tj(x,y) 



= mm 

0<k<2 



v{x,y) 



Following a discussion of Delone and Faddeev in [T] , we call the roots pi , p\ , p'{ 
of the equation F(x, a) = the left roots of the form F, while the roots p2, p'2, P2 
of F(d, —y) are called the right roots of the form F. If t\ is a left root, then it 
is easily seen that t% — —ad/tx is a right root of F. Two such roots of F will be 
called corresponding roots and we will assume that p\ and P2 , p[ and p' 2 , p'{ and 
p'2 correspond in pairs. 

The following lemma is a statement of Lagrange's method for solution of cubic 
equations by means of the resolvent adapted to the case of binary cubic forms. 

Lemma 3.1. For the cubic form F(x,y) the following identity holds 

1 



F(x,y) = 



where 



£ = £ix + &y, 
7] = 771 .t + n 2 y, 
£1 = pi + up[ + u 2 p", 
r)i = pi +uj 2 p' 1 +up", 

£,2= P2+ Up' 2 + UJ 2 P2, 
= P2 + ^ V2 + w /°2 

and cu = e~%~ . 

Proof. One can find the complete proof of Lemma 13.11 in [1] . 
We continue with the following definitions of p, q and itj : 

v + £ . V=i(v - 



□ 



p = 



V2 



q = 



V2 



(5) 



"1 



V6 



"2 



D-V«( « 
V6 



P 

Since 77 and £ are linear functions of x and y, so are p, q and it*. The reason for 
our interest in the new functions p(x, y), q(x, y) and Ui(x, y), despite their apparent 
complication, is that they explain the relation between the method of Evertse [5] 
and the method of Okazaki [15j for finding an upper bound for the number of 
integral solutions of (JT]) . In other words, these functions allow us to recast the 
resolvent forms £ and 77 in a geometric setting. 
By Lemma 13. 11 we have 



V6 



2^3 



2-v/3 



-[{Pi- P '{)x-{p2-p'l)y]. 
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g 

V6 



(p'l p'i> + (p' 2 - &)v 



We also have 

uj 2 - uj = cos(4tt/3) + v / ^Tsin(47r/3) - (cos(2tt/3) + v /Z Tsin(27r/3)) = ^=3. 
so we get 

(6) 
Further 

P 

V2 

(2pi + w[p^ + p'l) + u?(p\ + p'{)) x + (2p 2 + uj(p' 2 + p£) + oj 2 (p' 2 + gg)) y 

2 

Since w is a primitive third root of unity, u> + uj 2 = — 1 . Hence 

P _ 2 (Pl x + P2V) ~ {Pi + ~ (p' 2 + ^2)2/ 

v/2 2 



(7) 



Ml 


= D~ 


1/6 (P1 


-/>?)(*- 


P[y/a), 


"2 


= D~ 


' 1/6 (pi 


- Pi)(% - 


PxV/o), 


«3 


= D~ 






Piy/a), 



Substituting —ad/pi, —ad/p[ and —ad/p'{ for /?2> P2 an d P2 respectively, and noting 
that p\p\p'[ = —a 2 d, we obtain the following identities: 

(8) 
(9) 
(10) 

where p%, p\ and p" are left roots of F. Here we note that if we start with another 
choice of resolvent forms, only the order of Ui changes. In other words, all three 
resolvent forms can be indexed so that 

_ v^Tfa - &) 

*~ V2 

Let us assume that F is monic, as we may. Therefore 

(x - pxy)(x - ft x y){x - p"y) = F(x,y), 

If the pair (xo,yo) is a solution to ([I]), we conclude that [xq — piyo), (xq ~ p[yo) 
and (xo — PiVo) ar e units in Q(p\). Moreover, 

Ul u 2 u 3 = D- 1 / 2 { Pl - p'()(p[ - pM - p'{)F(x,y) = ±F(x,y). 

Suppose that (x, y) is a solution to F(x,y) — 1. Since 



log|ui| - log |w 2 ] = log 



and I }, v I is a unit, we can write 



Pi - Pi 



Pi - Pi 



log 



* - PiV 



(11) 



log - log |ua| = log Ai + mlogA 2 + nlogA 3 , 



where Ai = | p _^ | , A2 and A3 are fundamental units in the ring of integers of Q(pi) 
(when Z?i? > 0, the number held Q(pi) is real and has a ring of integer generated 
by a pair of fundamental units). 

Let us fix a resolvent forms (£i,?7i) and corresponding and g,. We get 
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By identities in (J3j) and Lemma 1370 \m\ — |&| and |£i(x, y)\ = y / H(x,y). Hence, 

Suppose that (x, y) is a solution to {!]) and related to resolvent form (£i,J?i) 
Since 



|l-£l = 



— = mm 

£i fc=l,2,3 



1 _ 2* 



we conclude that 



On the other hand, 



h 3 -£ 3 l 3^ 



2^ 2^2 

where the last equality comes from the equation j3}. 
If the solution (x, y) is related to (£,*, 77^) , then 



(12) 

So we have 



|« 3 | = \D- 1/6 ^Qi\ < 1. 



log|it 3 | < 0. 

The identity 

|wi(x,y)u 2 (x,y)u 3 (x,y)| = 1 
holds when (x,y) is a pair of solution to \F(x, y)\ = 1. Therefore, 

log I tti| + log|u 2 | +log|u 3 | = 

and 

log |miu 2 | > 0. 

4. Geometric Gap Principles 

We will study the geometric properties of the functions Ui defined in section 2, 
by considering the well-known geometric properties of the unit group U of Q(pi), 
where p\ is a root of F(x, 1) = 0. 

Since we assumed that F has positive discriminant, the algebraic number field 
Q(pi) is real and has two fundamental units, say A 2 and A3. By Dirichlet's unit 
theorem, we have a sequence of mappings 

U 1 — ► V — {0} c M 3 

and 

log:F-{0}^ A-{0}, 

where A is a 2-dimensional lattice, U 1 — > — {0} is the obvious restriction of the 
embedding of K in M 3 , and log is defined as follows: 
For (x 1 ,x 2 ,x 3 ) G V"- {0}, 

log(xi,x 2 ,x 3 ) = (log|xi|,log|x 2 |,log|x 3 |). 
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We define r to be the embedding from the unit group U to the lattice A: 

t : U i — ► A - {0}. 
By identities ©, © and ([ID)), the vector 

tt = (log |«i | , log |« 3 1 , log |tt 3 1) 

can be considered as 

(13) v + (log \x - p' t y\, log | a; - p"y|, log \x - p x y\), 

where 

= (log \D~ 1/6 ( Pl - p'{)\ , log II?- 1 / 6 (p' x - pi)|,log|£>~ 1/6 (^ - 

We have assumed that F(x,y) is monic, so we can suppose that (1,0) is a pair 
of integer solutions to F(x, y) = 1. Note that the vector v in (fT3| is a permutation 
of the vector u(l, 0). 

If (x, y) is a solution to |F(x,y)| = 1, then 

u e v + A = Ai. 

Note that Vol(A) = Vol(Ai), where Vol(A) is the volume of fundamental paral- 
lelepiped of lattice A. Since u belongs to a 2-dimensional lattice, we can find a 
2-dimensional representation for it, say (t,s). Specifically, let (x,y) be a solution 
to F(x, y) — 1 and define functions t and s of x and y as follows 



t = 



-V6 



log |w 3 | 



log I "1 1 - log I "2 I 



Then we have 
(14) 



log = s/V2 + t/V6 
log |u 2 | = -s/V2 + t/V6 
loglual = -2t/V6. 
Therefore, it can be easily verified that 

u= (log |ui|, log^l) logKI) = sa + t/3, 
where a = ^=(1,-1,0) and j3 = ^=(1,1,-2) are two orthonormal vectors in R 3 . 

Hence, we can write u = (t, s) and ||u|| = Vs 2 + 1 2 , where || || is the L 2 norm. By 
(HH) we get 

( log ^1 ,log ^ ,log ^ ) = Visa' + V3tP', 

\ U 2 U 3 Ul J 

where 

1 -1 -1 

\/3 V ' VH %/6 
Since a' and /?' are orthonormal vectors in R 3 , we get 



a' 



and 



(15) 



flog 




,log 


u 2 


,log 


^3 


) 




u 2 


U3 







V3\/.s 2 + t 2 = V3|| 



Remark. Since log \ us\ < 0, the function t is a positive- valued function. 
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Lemma 4.1. If s > then log|ui| > log|ii2| and 

2sinh(s/v2) = exp(— Vet/2), 

and if s < then log|ui| < log |ti2 1 and 

2 sinh(— s/v 2) = exp 

Proof. From (JSJ), the definition of Ui, we have u\ +112 +W3 = 0. Assume that s > 0. 
From (|f4p . since s and t are both nonnegative, we get \u\\ > |«2| and |wi| > \v,3\. 
Therefore, 

e s/V2+t/V6 _ e - s /V2+t/Ve _ e -2t/Ve = | Ul | _ |„ 2 | _ | U3 | = , 

and 

Noting that e s /^ — e ^ s /V2 — 2 sinh(s/ v / 2), will complete the proof. One can give 
a similar proof for negative s . □ 

Let us define 

,(,) := VSstoh- (^Et^I2l"j . 

Then s = ±g(t). 

In the following theorem, we summarize the properties of function g, which will 
be used later. 

Theorem 4.2. Let g(t) = v^sinh -1 f ^st^MM \ . We have: 

(i) g is decreasing . 

(ii) For any t > 0, 

\s\ = g{t) < e-^/ 2 /V2. 

(iii) The function g(t)e at is decreasing when a < 

Proof. (i) Since 

sinh( ff /V2) = exp(-V6i/2)/2, 
we have the following implicit differentiation: 

Jcosh( 3 /V2) = Z^e*p(-V6*/2). 
Since cosh(g/\/2) and exp(— y&t/2)/2) are both positive, 

^<0. 

(ii) Define the function 

f{x) = v / 2sinh(a;/\/2) - x. 

The first derivative test shows that / is an increasing function and for 
positive x, f(x) > /(0) = 0. So 

\/2smh(a;/v / 2) > ac, 
when x > 0. Put a; = |s| to get 

\s\ =g(t) < V2smh(\s\/V2). 
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(iii) Set 

A{t)=g{t)e a \ 

then 

A'(t)=e at (g'(t)+ag(t)). 
For a < 0, A' < since g' < 0. For positive a, by part (i) and (ii), we have 

— a \ 



A'(t) < expat- V6t/2 



2cosh(g/V2) \/2 



Since g is a decreasing and positive- valued function, cosh(<?(t)/V2) is a 
decreasing function of t. So we have 



cos 



h(g{t)/V2) < cosh(<?(0)/V2). 



An easy way to evaluate cosh(<?(0)/v2) is to recall that 
smh(g(0)/y/2) = exp(0)/2 = 1/2. Therefore, 



cosh( 5 (0)/V2) = Jl + i = V5/2. 

We conclude that A' is negative if 

-V3/V5 + a/V2 < 0. 

This means 

a < —=. 
y/5 



□ 



Lemma 4.3. Let (x,y) and (x',y') be two distinct solutions to equation (Qp, related 
to (77, £). Put p = p(x, y), p' = p(x',y'), q = q(x, y), q' = q(x' , y'). We have 

\pq' -p'q\>VSD. 

Proof. By definition 

V + £ V=ltf - £') rf + e V^T(r; - 



\pq' -p'q\ 



V2 V2 V2 V2 



Since ^(x,y)r](x,y) = H(x,y) is a quadratic form of discriminant —3-D , it follows 
that 



r£ - i/C = ±V^3D(xy' - x'y). 

Since (x,y) and (x',y') are distinct solutions to F(x,y), xy' — x'y is a nonzero 
integer. □ 

Lemma 4.4. Le£ (x, y) and (x',y') be two distinct solutions to equation (Qp, related 
to {£,T)). Assume that t(x' ,y') > t(x,y). Then we have 

t(x\ y') > 2t(x, y ) + ^-\ogD- V6 log ( 2 + ■ 
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Proof. Put 



and 

First we show that 

By the triangle inequality we have: 
\p\ = -±=\V2p\ < 



p 


= p{x,y), 


P 


= p{x',y'), 


q 


= q(x,y), 


q 


= q(x',y'), 


s 


= s(x,y), 


s' 


= s{x' ,y'), 


t 


= t{x,y), 


t' 


= t(x',y') 


Ui = 


= Ui(x,y), 


u 


= Ui{x ,y ) 


\p\ 


< V2D 1 / 6 


e t/VS 


cosh(s/\/2). 



L (|(p/V2) + (q/V6)\ + \(g/V6) (p/V2) 



+ 



V2 

1 \(g/V6) + ( P /V2)\^ 
V2\(q/V6)-( P /V2)\^ 
1 \(q/V6)-(p/V2)\^ 



1/2 



(q/V6) + (p/V2) (q/y/6) - ( P /V2) 



1/2 



(q/y/6) + (p/V2) (q/V6) - (p/V2) 



1/2 



1/2 



V2\{q/V&) + {p/^/2)\^ 
By (Jill) and ((231), we have 

(q/V6) + {p/V2)\ 1/2 \(q/V6) - (p/V2)\ 1/2 = D^u^u^ 2 = D^expit/VG). 

Equations fT4")) and (f2"5)) also give us the following identities: 

( \(g/V6) + (p/V2)\^ \(q/V6)-(p/V2)\^ 
\\(q/V6) ~ (p/^V 2 \(q/V6) + (P/V2)\^\ 

(lo 8 |ui|-log|^ 2 l) (log -log |»i2 I) 

= e 2 -f e 2 



2cosh(^=) 



and 



\q\ = (V6/2)^ 1 /6 e -2 t /76^ 

Using Lemma POl we get 

D i/ e < e (t'-m/Ve cosh n s/ y^ +e (t-^')/Ve cosh n s y^\ 

One can express the above equation in terms of sinh instead of cosh by substituting 
cosh(|s|/V2) with 

smh ^ifT) + e -M/A 

Now we use the assumption that t' > t and the fact that er\ s \l^- < \ gy Lemma 
14.11 we get 

/ \ / „-Vet/2\ 

D 1/6 < e^'- 2 *)/^ (i + e -3(t'-t)/Vej i + e 
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Note that by Theorem (14. 2| . t > log(2)/v / 6, whereby taking the logarithm of both 
sides of the above equality, yields 

t'-2t> V6/6\og{D) - V61og ^(l + e -3(*'-*)/^ (\ 

Therefore, 



2V2 



t' - 2t > V6/61og(L>) - V61og ^2 + i 



□ 



Lemma 4.5. Suppose that (Qp /ias i/iree distinct solutions related to (£,77). TTien 
i/iree distinct corresponding points (t, s) , (t',s') and (t",s") form a triangle. 

Proof. Suppose (t, s) , (f , s') and (t", s") are collinear and t < t' < t" . Then 

s' -s s" - s' 



t' -t t" - f 

Assume, without loss of generality, s' > 0. Since g(t) = \s\ is a decreasing function 
of t, we have 

s" - s' < 
and consequently, s' — s < 0. Therefore, 

s > 0. 

Since we assumed (t, s) , (t' , s') and (t", s") to be collinear, 

s"-s s'-s 



t" -t t'-t 

By Lemma B~4l and since \s\ > \s'\ > \s"\ > 0, we get 

s'-s -s' -2s' s"-s' 

—t < — < — < — < °- 
This contradiction shows that (t,s) , (t',s') and (t",s") are not collinear (note 
that any vertical or horizontal line intersects the graph of g and g' at most in two 
points). □ 

Suppose that |T]) has three distinct solutions related to (£, rf) and A is the area of 
the triangle formed by three distinct corresponding points (t, s), (if , s') and (t" , s"). 
Then vectors (t — t',s — s') and (f — t" , s — s") generate a sub-lattice of Ai with the 
volume of fundamental parallelepiped equal to 2A. Therefore, 

2 A > Vol(Ai). 

Now let us estimate 2A, the area of rectangle which has (t, s) , (t', s') and (i", s") as 
three of its edges. Recall that s(x,y) = ±g(t(x,y)) and g is a decreasing function. 
Suppose that t<t'< t" . Then g(t") < g(t') < g{t) and we have 

2A<(t"-t)(g(t)+g(t')) = (t"-t)(\s\ + \s'\). 

Part (iii) of Theorem 14.21 shows that 



-V6(t' -t) 

Is'l < Isle ^ 



Therefore, 



/ -V6(t'-t) \ 

Vol(A) = Vol(Ai) < (t" - t)\s\{ 1 + e J 
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Using Theorem 14.21 again, we get the following gap principle of this paper which is 
essentially Theorem 5.5 of [15] : 

Theorem 4.6. Suppose that F(x,y) has three distinct solutions (x,y), (x',y') and 
(x",y"), all related to (£,tj). Assume that 

t = t(x,y)<t'=t(x',y')<t"=t(x",y"), 
where t is the function defined in the begining of this section. We have 

„ > V2Vol(A)exp(VEt/2) 
~ l + exp(-\/6(t'-t)/v / 5)' 
where Vol(A) is the volume of fundamental parallelepiped of lattice A. 

5. Linear Forms In Logarithms 

We have seen that \/2s = log — log |m 2 | = log Ai -fmlog A2 + nlogA3. Where 
s is a function of (x, y) defined in Section 3 and m are also functions of (x, y) defined 
in Section [2] By Lemma |4~2"1 we have 

log(V2|s|) < -(V6/2)i. 

Here, we will use a well-known lower bound for linear forms in logarithms of alge- 
braic numbers, to find an upper bound for log(v / 2|s|). 

Theorem 5.1 (Matveev). Suppose that K is a real algebraic number field of degree 
d. We are given numbers ai, . . . a n € K* with absolute logarithmic heights h(aj). 
Let logai, . . . , loga„ be arbitrary fixed non-zero values of the logarithms. Suppose 
that 

Aj > m&x{dh(aj), | loga^l}, 1 < j < n. 

Now consider the linear form 

L = 61 logai + . . . + b n log a„, 

with bi, . . . , b n € Z and with the parameter B — max{l, max{bjAj/A n : 1 < j < 
n}} . Put 

n = A 1 ...A n , 

1 fi 1 

C(n) = —e n (2n + 2)(n + 2)Un + 4) n+1 (-en), 
n\ 2 

C = log(e 4 - 4n+ V- 5 d 2 log(en)), 
W = log(1.5eSdlog(ed)). 

If b n 7^ , then 

log I L I > ~C{n)C Wod 2 n. 

Here, we recall the definition of absolute logarithmic height from [HI [TU] . Let 
Q(p) cr be the embcddings of the real number field Q(p) in R, 1 < a < 3, where p is 
a root of F(x, 1) = 0. We respectively have 3 Archimedean valuations of Q(p): 

|a| a = |aW|, l<a<3. 

We enumerate simple ideals of Q(p) by indices a > 3 and define non- Archimedean 
valuation of Q(p) by the formulas 

\a\ a = (Norm p) _fe , 

where 

k = ordp(a), p = p a , a > d, 
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for any a £ Q* (p) . Then we have the product formula : 

oo 

JJ \a\cr = 1, a € 



1 

Note that \a\ a ^ 1 for only finitely many a . We define the absolute logarithmic 
height of a as 

_j oo 

h(a) = -J2\l°E\*\<r\- 

We will apply Matveev's lower bound to 

log — log |m 2 | = log Ai + mi log A 2 + ri\ log A 3 . 

Suppose that 

\\u(x Q ,y )\\ = min ||w(a;,y)|| 

and 



log|ui(x ,?/o)| ~ log|u2(a;o,yo)| = log 



P-P 



P' 



a log Ai + 6 log A2 



then for any solution (x,y), we can write 

log|ui(x,y)| - log\u 2 {x,y)\ = logA + mlogAi + nlogA 2 , 
where m = mj — a, n — n\ — a and 



A 



P-P" 



P ~ P 

Since A 2 and A3 arc the fundamental units of the ring of integers of Q_(p), Xi, A 2 
and A3 are multiplicatively dependent if and only if Ai is a unit. If Ai is a unit then 
we can write log \ui \ — log |u 2 | as a linear form in two logarithms. Since Theorem 
15.11 gives a better lower bound for linear forms in two logarithms, we can assume 
that Ai, A 2 and A3 are multiplicatively independent and log \u±\ —log |u 2 | is a linear 
form in three logarithms. 

First, suppose that A is a unit in the number field. We have 

MA) = \ (|log|A|| + |log|A'|| + |log |A"||) = \ |r(A)| x , 

where A' and A" are the conjugates of A , r is the embedding of units to the lattice 
A and | |i is the L\ norm on R 3 . So we have 

MA) = i|r(A)| 1 <^||r(A)||, 

where || || is the L 2 norm on R 3 . So when A is a unit 
(16) max{3/i(A),|log(|A|)|} < ||r(A)||, 



since |log(A)| < ^/log 2 |A| + log 2 | A'| + log 2 |A"| = ||r(A)||. 
In the identity 

log |«i I -log|u 2 | = log A + mi log A 2 +?iilogA 3 , 

A 2 and A3 are fundamental units of Q(p). Therefore, in Theorem 15. 1[ Ai can be 
taken equal to |r(Ai)| 2 , for % = 2, 3. 

Bases h\ and £> 2 of lattice A are called reduced if the following conditions are 
satisfied : 



CUBIC THUE EQUATIONS 



15 



(i) bi < \\v\\ for every vector iJeA - {0}; 

(ii) &2 < \\v\\ for every vector v G A — Zb%. 

Remark. Although the definitions of reduced basis for lattices and reduced forms 
are somehow related, one should note that we define them separately and they are 
not to be confused. 

It is a fact that we can always choose a pair of reduced basis for a two dimensional 
lattice. So we choose the fundamental units A2 and A3 such that the basis r(A2) 
and r(Aa) are reduced basis for A. When 61 and 62 are the reduced basis of A, 
since 61,62 < 61 ± 62, we conclude that the angle between vectors 61 and 62 must 
be between tt/3 and 27r/3. Therefore, A2 and A3 can be chosen so that 



||r(A 2 )|| ||r(A 3 )|| < ^Vol(A). 



Hence, in our case, 



By (fT5)), we have 



A 2 A 3 < -=Vol(A). 
v3 



n \ u i\ , F2 , "3 

(log 1 r, l0g— , log 1 r 

"2 \U3\ \Ui\ 



V3\\u\\. 



The well-known inequality ^±§±2 < [2!+^+^] 1/2 s h 0W s that 

(17) < V3\\v\\ 

for every vector »el 3 . Therefore, 

|log \(p - p"){x - p'y)\\ + |log \(p' -p)(x- p"y)\\ + |log \(p' - p"){x - P y)\\ 
< 3 ||tt(x 0) yo)|| • 
Now, we note that 

(p-p")( X -p'y) 



E 



cr>3 



log 



(p> - p){x- p"y) 



^2\log\(p'-p)(x-p"y)\J. 



cr>3 



cr>3 



We know that Archimedean valuations of \(p — p"){x — p'y)\ are itself, 

\(p' — p"){x — py)\ and \p' — p(x — p'y)\. So by the product formula, since (x, y) is 

a solution to ([T|), the product of all non- Archimedean valuations equals D~ x l 2 . 

Therefore, 



J2log\(p-p")(x-p'y)\ c 



CT>3 



and similarly 



^log|(p'-p)(a;-/y)| c 



(T>3 



= \logD, 



= - 2 lo g D. 



Since (p' — p)(x — p"y) and (p — p"){ x ~ p'y) are algebraic integers, we get 



(18) 

This gives an estimate for A\ 



h(\i) < -(3\\u{x ,yo)\\ + log D) 
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Let B\ = BA3, where B is as in theorem 15. II Then 

B x = max{bjAj, : 1 < j < 3}. 

Since 



log 



we can write 



log 



,log 



= log 
,log 



{p-p"){x-piy) 



(p' -p)(x- p"y) 



m\i + n\ 2 , 



\{p-p"){x-p'y)\ Up>-p)( x - p "y)\ \{p/-pf'){x-fyu)\ 

g \{p>-p){ X -p»y)y g \{p> - p»){x - py)V g \(p-p»)( X -f/y)\ 

+ m\i + n\ 2 , 

where A.; = t(A^) , for i = 2,3. Since A2 and A3 have been chosen so that A2 and 
A3 form a reduced basis for the lattice A, we get 

m|A 2 |i ,n|A 3 |i < 



log 
log 



\(p- 


-p"){x- 


p'y)\ 


,k>£ 




- p)(x- 


\{p' 


-p){ 


x — 


P"y)\ 


5 l(p' 


- P")(x 


Ui 


,log 


U2 


,log 


U3 


) 






U 2 


«3 


Ui 




1 





plog- 



Therefore, by (T5D and ([T7J) 

(19) mpf 2 |i,n|Al,|i <3(||u|| + \\u(x ,y )}\) . 

Theorem 5.2. Let F be a cubic binary equation with positive discriminant. For 
all pairs of solution (a, b) to the equation (QP, except possibly one of them, we have 

D < Me 2VEt , 

where D is the discriminant of F{x,y) and t = t{a,b), for the function t defined 
in Section^ Moreover, when t > 5 

1 



D < -e 



2V6t 



Proof. By © 

By ©, 
and 



\h\ = N = 



,p_+iq p-ig , = p +g 



\ P \ = \^( Ul - u 2 )\d^ < ^V/ 6 (M + |u 2 |). 

Therefore, by (JTU) 

(20) |ff I < -D 1 / 3 ( e 2*/y6 (e 2 S /y2 + e -2s/V2 + 2 )/2 + -e- 4 */ 6 



Therefore, by Lemma |2~T| for all solutions (a, b) to (p}, except at most one of them 
(21) ^ 2 V3 < i^ 1/3 ^ e 2t /V6 (e 2,/y2 + e -2s/V2 + 2 )/2 + V 4t / 6 ^) . 
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Since t > 0, part (ii) of Lemma Q] says that \s\ < e _v/ ^'/ 2 /\/2- Hence, 

£,1/8 < 2e 2t/ ^, 

which proves the theorem for general t. When t > 5, we note that f e" 4 */ 6 < 0.054 
and |s| < 0.0016 by Theorem [42] . " □ 



Since u — (t, s) and ||u|| = \/t 2 + s 2 , from Theorem 14. 2[ we deduce that \u\ 2 is an 
increasing function of t. So we can assume that Theorem 15.21 is satisfied for all 
solutions, except possibly (xo,yo), where 

ll"(a;o,yo)|| = , min \\u(x, y)\\, 

S is the set of all solutions to (QJ and u = (log log \v,2\, log \u^\). 

Suppose that three distinct solutions (x,y), (x',y') and (x" ,y") of (P) are related 

to (£,77) and t" = t(x",y") > tf = t(x',y') > t = t(x,y). First, we recall that 



u = (t, s) and \u\ 2 = \Jt 2 + s 2 . By Theorem l4~2l if wc take t > 5 (Theorem IS~2 
enables us to assume that t is large), we get 



u 2 \/t 2 + e-^6 t /2 < 1.00161 
Therefore, by (fT5]l and theorem 15.21 we can take 

A 1 = l.Q016(|t + V6t). 

Inequality (31]) suggests the value 3(1.0016)(i" + 1) for B 1 . But by theorem IP1 
for large discriminant D, t" > At . So we take 

(22) B x = 3(1.0016) (*" + t"/A). 

So Matveev 's lower bound gives us: 

log|i| > -1.5036 x 10 11 AiA 2 A 3 log(25.6708(3.0048)(1.25)t7A 3 ) , 

where L = log \u\\ — log \u 2 \. On the other hand, by Theorem l4.21 we have 

log|£| = logv^s" | < -VQt"/2. 

We conclude that 

t" < 1.2276 x 10 11 A 1 A 2 A 3 \og(96.27blt"/A 3 ), 



or 



Therefore, 



96.2751*7* < L1892xl0 i3 AiA2 . 



log(96.2751i"/A 3 ) 



log(96.2751i"A4 3 ) < log(1.1892 x lO 13 **). 

e — 1 

Recalling that A 2 A 3 < 2/ v / 3Vol(A), we obtain the following upper bound for t": 

(23) t" < —^—1.2276 x 10 11 (4= I VolfAUi log(1.1892 x 10 13 AiA 2 ). 
e-1 V v3/ 
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6. Proof Of The Main Results 

Let 

\\u(x ,y )\\ = min {Xty)eS \\u(x,y)\\. 

Suppose that (x,y), (x',y') and (x",y"), with none of them equal to (xo,yo): are 
three distinct solutions to (JTJ) , and related to a fixed choice of resolvent form. Let 
t = t(x,y) < f = t(x',y') < t" = t(x",y"). By dHJ and Theorem|4~H we get 

' l -1.2276 x 10 1 M 1 (-^=)Vol(A)log(1.1892 x lO 13 ^/^) 



e-1 ' V V3 
> y / 2Vol(A)exp(V6V2) 
" 1 + exp(-\/6(i' - t) /VS) ' 

where A 1 = (3/2 + V6)(1.006)i and A 2 = ||r(A 2 )||. Without loss of generality, we 
can assume that ||r(A2)|| < ||r(A3)||. Therefore, 

^|r(A 2 )|I<^||T(A 2 )||||T(A 3 )||<Vol(A). 

We have 

log|wi(a;,y)| - log|tt 2 (a;,y)| = logAi + m' log A 2 + n' log A 3 , 

where m! and n' are integers. Since (x, y) ^ (xq, yo), at least one of m' or n' is a 
nonzero integer. So by (fT9j) , we have ||t(A 2 )|| < 6.01i . Using Theorem 14. 4i we get 

4.8484 x 10 n (-^=)Vol(A)ilog(2.8184 x 10 14 t 2 ) 



1-e ' V V3 

> %/2Vol(A)ea;p(V6t/2) 

~ 1 + cxp(-v / 6t/v / 5) 

Therefore, t < 27.91 and by equation (f2Tj) , 13 < 5.31 x 10 59 ; i.e. we have proven 
that there are at most 2 pairs of solutions (x, y) ^ (1, 0) and (x' ,y') ^ (1, 0) 
related to a resolvent form (£, rj) , when D > 5.31 x 10 59 . 
If we suppose that D < 5.31 x 10 59 , by (|18p. we can take 

A 1 = ^+ilog(5.31xl0 59 ). 

By substituting this new value of A\ in (f2"3")l . we get 

t < 27.5321, 

and therefore, D < 1.4 x 10 57 . Since we have three pairs of resolvent forms, 
Theorem (11.11) is proved. 

As we mentioned in the remark after the proof of Lemma l2.ll the solution (1, 0) 
needs to be treated separately, only if F is equivalent to a monic reduced form. 
Otherwise, (xo,yo) ^ (1,0) and Lemma I2TT1 and therefore Lemma IS~2l will hold for 
all solutions without any exception. By the analytic class number formula and 
Louboutin's upper bound (which can be found in [3]) : 

Vol (A) < ^^log 2 ^, 
8 

we have that 

A 2 < i^ 1/4 logZ3. 
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By Theorem E21 

A 2 <I e v^/> g I + 2Vtt). 

Now, having appropriate values of A\ and A2 in hand, we solve inequality (|24[) to 
get t < 28.38 and consequently by (f2~Tj) . fl< 9 x 10 58 ; i.e. we have proven that 
there are at most 2 pairs of solutions (x, y) and (x',y') related to a resolvent form 
(£, 7]) , when D > 9 x 10 58 . Therefore, we get Theorem 1 1.2 1 

In [2] , it is proved that if D > 2400, related to a fixed pair of resolvent form, there 
are at most 3 different pairs of solutions {x,y) to |T]) with H(x,y) > |\/3Z), where 
H is the Hessian of F. This together with lemma |2~T1 leads to the main theorem of 
[2], that is, the equation F(x,y) = 1 has at most 10 solutions in integer x and y . 
For < D < 2400, equation F(x, y) = 1 is completely solved for representatives of 
every equivalent class of binary cubic forms. These computations show that the 
equation |T]) with discriminant < D < 10 6 has at most 9 solutions in integers x 
and y. The complete result of these computations are tabulated in section 9 of [2J. 
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